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Abstract 

Ge, Rusjan, and Zweifel (J. Stat. Phys. 59, 1265 (1990)) introduced a binary 
tree which represents all the periodic windows in the chaotic regime of iterated 
one-dimensional unimodal maps. We consider the scaling behavior in a modified 
tree which takes into account the self-similarity of the window structure. A non- 
universal geometric convergence of the associated superstable parameter values 
towards a Misiurewicz point is observed for almost all binary sequences with 
periodic tails. There are an infinite number of exceptional sequences, however, 
which lead to superexponential scaling. The origin of such sequences is explained. 



1. INTRODUCTION 



Iterated one-dimensional unimodal maps [1] have been the subject of hundreds 
of papers over the last couple of decades. Not only are they among the 
simplest dynamical systems exhibiting chaotic behavior but also very important 
as prototypes of dissipative systems. In spite of breakthroughs such as symbolic 
dynamics [1,2,3], ergodic behavior [1,4,5], and the transition to chaos [6] in the 
theory of unimodal maps, many problems remain related to the scaling behavior 
of such maps within the so called chaotic regime. This region is defined as the 
parameter interval between the first period-doubling accumulation point and the 
final crisis point beyond which no periodic or chaotic attractors can be found 
within the unimodality interval of the phase space. 

Rigorous mathematical proofs [5] establish that the parameter values correspond- 
ing to an absolutely continuous invariant ergodic measure form a set with a positive 
Lebesgue measure. These "chaotic" parameter values are found in between the 
infinite number of windows with stable periodic attractors. The "periodic" win- 
dows are expected to be dense on the parameter axis. Although each window has 
a finite lenght there remains a great deal of space for chaotic parameter values: 
Near the accumulation point of a period-doubling cascade, the relative fraction of 
the aperiodic solutions is given by the universal number 0.892... [7]. Even consid- 
ering the whole chaotic region, the probability of finding an aperiodic solution is 
approximately 9/10 for a typical map [8]. 

Since periodic windows are ubiquitous along the parameter axis, various infinite 
sequences of them are a natural tool when investigating scaling properties of 
unimodal maps. The scaling behaviors of period-doubling [6] and more general 
multifurcation sequences [9], period-adding sequences approaching a crisis [10,11] 
and tangent bifurcation points [12] have been determined. Shibayama [13] 
extended the analysis to the so called Fibonacci sequences whose scaling is 
superexponential both on the parameter axis and in the phase space. An exact 
universal form of that type of scaling was found by Ketoja and Piirila [14] using a 
renormalization argument. Later Lyubich and Milnor [15] derived rigorous results 
both for the scaling behavior and the dynamics at the accumulation points of such 
sequences. 

In this paper a multitude of new examples of both non-universal geometric and 
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universal superexponential scaling are reported and the corresponding sequences 
related to a binary tree of periodic windows introduced by Ge, Rusjan, and Zweifel 
(GRZ) [16]. Originally the tree was defined so that each window included the 
period-doubling tail in addition to the stable parameter interval of a periodic 
solution. In our modified tree, each window is extended up to the corresponding 
interior crisis point. We "sum" not only over the period-doubling tail but over 
all the multifurcation sequences. In this way the self-similarity of the periodic- 
window structure [17] can be naturally taken into account. The structure within 
each window is essentially a small copy of the entire structure along the parameter 
axis. We concentrate on those periodic-window sequences whose binary codes have 
periodic tails. Such sequences usually lead to arithmetic growth of the period and 
to non-universal geometric scaling. The windows accumulate at a Misiurewicz 
point [1,4] at which the dynamics is completely chaotic. In addition, we find an 
infinite number of exceptional cases in which the period increases geometrically and 
the scaling is superexponential generalizing earlier results on Fibonacci sequences. 
Finally, preliminary numerical conclusions on aperiodic binary codes are reported. 
It is conjectured that the accumulation point of such a sequence always corresponds 
to a chaotic attr actor. 

2. BINARY TREE 

Consider a one-parameter family fn{x) of differentiable unimodal maps from a 
real interval / to itself. is assumed to have a quadratic maximum at a; = c 
so that the map is monotonically increasing for a; G /, x < c and monotonically 
decreasing on the other side of the critical point c. An orbit obtained by iterating 
starting from c can be symbolically represented in terms of the kneading 
sequence aia2... where ai = R \i fjj_{c) > c and ai = L if /^(c) < c. The case 
in which the orbit returns back to the critical point after i iterations, /^(c) = c, is 
indicated by cutting the kneading sequence after i — 1 symbols so that the sequence 
becomes finite. Finite symbol sequences therefore correspond to superstable 
periodic orbits. We are interested in the admissible kneading sequences at some 
parameter values of the unimodal map. Metropolis, Stein, and Stein (MSS) [2] 
discovered a simple rule by which all admissible symbol sequences, the so called 
MSS sequences, can be generated and arranged on the parameter axis. Originally 
the rule was developed with one-parameter families of one-dimensional maps in 
mind, the type where the parameter appears as a multiplicative factor in the 
definition of the map. It holds, nevertheless, in a much larger class of unimodal 
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maps, e.g. the logistic map fij,{x) = ji — x^. 



It is instructive to take a look at the origin of the rule. To this end, consider 
two parameter values ni and 112 (a*i < A*2) which correspond to the infinite MSS 
sequences A and B. As in ref. [16], the beginning shared by the sequences is 
denoted by Af\B. In other words, //i corresponds to the symbolic orbit {Af\B)an--- 
and /i2 to the orbit {A/\B)hn---, where the symbols a„ and differ. The parameter 
dependence is assumed such that the MSS sequences in the interval (/ii, ^12) always 
begin with A /\ B. By continuity there has to be a parameter value /^s G ^2) 
with the finite symbol sequence C = A/\B. The corresponding orbit is superstable 
with the period n. Within the interval {pi,H2), there are no other orbits with 
periods < n. Above /X3, but within the window for the stable period n, the MSS 
sequence has the form h{C) = CbnCbnCbn... (the nth symbol has to be bn — 
otherwise there would be at least two parameter values with the MSS sequence 
C). From h{C) and B a new periodic window in between and /i2 can now be 
constructed. Similarly, the infinite sequences A and a{C) = CanCanCan--- can be 
chosen to generate a periodic window in between /ii and ^13. The only way this 
recursive procedure of generating new periodic windows can get stuck at some 
level is two infinite MSS sequences' becoming equal. By repeating the procedure 
an infinite number of times, MSS sequences can be generated for aperiodic orbits 
as well. 

One needs two MSS sequences in order to get the procedure going. For a so-called 
full one-parameter family of unimodal maps [1] , such as the logistic map, one can 
set out with the symbol sequences for the superstable period-two cycle (R) and 
the final crisis point (RLLL...). In the latter, c is mapped onto an unstable fixed 
point after two iterations. 

The infinite sequences h{C) and a{C) above are called the harmonic and 
antiliarmonic extensions of C. These extensions can be expressed without knowing 
the infinite "parent" sequences A and B. Let us first write the extensions in the 
form 

h{C) = CaCaCa... 
a{C) — CaCaCa... 

where a denotes the "conjugate" of the symbol a; i.e. R = L and L = R. In the 
sequel we will refer to a or a as binding elements of an extension. Around the 
critical point c, there is a small box in which the nth iterate of the unimodal map 
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looks very similar to the first iterate. The nth iterate has either a maximum or a 
minimum at c. At the superstable parameter value with the MSS sequence C, the 
extremum touches the critical point. As the value of the parameter is increased 
within the stability interval, the extremum passes either below or above the critical 
point depending on whether the extremum is a maximum or a minimum. If C 
contains an even number of ii's (C even) the nth iterate has a maximum at c; 
otherwise {C odd) the extremum is a minimum [11]. In other words, a = i? in the 
former and ck = L in the latter case. 

The fact that the nth iterate restricted to a small box around the critical point 
becomes a unimodal map is responsible for the self-similarity of the MSS structure. 
One expects the same MSS periodic windows with the nth iterate as with the 
original map. Only the structure in the higher iterate is observed in a much smaller 
parameter interval than for the original map. The window with the "reduced" MSS 
sequence ai...afe, which includes only every nth iterate (these actually land inside 
the small box), corresponds to the full MSS sequence [17] 

„ . ^ — / CaiCa2---CakC, if C is even 

G * {ai...ak) - j CaiCa2...CakC, if C is odd 

for the original map. By this composition law it is easy to write down the MSS 
sequence at the endpoint of the parameter interval which contains the self-similar 
copy of the whole periodic window structure: C * (RLLL...) = CaCaCa... = 
Caa{C). This endpoint corresponds to an internal crisis where the orbit of the 
critical point lands on an unstable period n after 2n iterations. We call Caa{C) 
the crisis extension and denote it by e(C). 

The superstable period-two cycle with the MSS code R, preceded by a stable 
fixed point on the parameter axis, belongs to the primary period-doubling cascade 
which ends at the transition to chaos. By self-similarity, every periodic window 
is followed up by a similar cascade. GRZ [16] modify the definition of a periodic 
window including in it the corresponding period-doubling tail. In the recursive 
construction of the periodic windows one considers, instead of h{C), the MSS 
code at the period-doubling accumulation point, h{h{h{...h{C)...))), where the 
"cut" harmonic extension h{C) = CaC is successively applied an infinite number 
of times. In this way, one never generates windows within a period-doubling 
cascade and two neighboring infinite sequences never become equal. Therefore, it 
is possible to generate an infinite binary tree of periodic windows. GRZ want to 
apply Feigenbaum's [18] general ideas on the renormalization of binary trees to 
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this particular case. We take a different point of view and exploit the binary tree 
as a tool of studying the overall scaling behavior within the chaotic region. 

We modify the GRZ tree so as to take into account the self-similarity of the 
window structure. Each window is extended up to the corresponding interior crisis 
point. The infinite binary tree is constructed taking advantage of the antiharmonic 
and crisis extensions in the following way: 

1) Begin with the infinite "left parent" e{R) and the infinite "right parent" 
RLLL.... The crisis extension e{R) = RLRRR... corresponds to the last band 
merging point (where a critical point is mapped onto an unstable fixed point after 
three iterations) and the sequence RLLL... to the final crisis point. All the windows 
of the tree lie in between these two points (i.e., in the region of one chaotic band 

[I])- 

2) From two infinite "parent" sequences A and B form a finite "daughter" sequence 
C by taking the shared beginning of A and B. The first such sequence RL is the 
"root" of the tree. 

3) Take A and a{C) as the infinite parents of a new left "branch" and e(C) and 
B as the infinite parents of a new right "branch". Attach the symbol to the left 
and the symbol 1 to the right branch. 

The beginning of the infinite binary tree is shown in Fig. 1. There is a one-to- 
one correspondence between the binary codes consisting of the symbols and 1 
and the MSS sequences. In the following, we use the symbol to express this 
correspondence. For example, 10^ — > RL^RLR. a'' means that the symbol (or a 
block of symbols) a is repeated k times. This convention is used for both the MSS 
sequences and the binary codes. 

3. TRANSFORMATION BETWEEN THE BINARY CODE AND 
THE MSS SEQUENCE 

In this section we develop " self-contained" recursive rules by which the transfor- 
mation — > can be carried out. These rules are just another variant of the MSS rule 
but they turn out to be the key to understanding how the the lenght of the MSS 
sequence increases down the binary tree. In the following, the ith. symbol of the 
MSS sequence A is denoted by {A}i and the string from the zth symbol up to the 
jth symbol by {A}1 [j < i implies an empty string) . This notation is particularly 
useful if A is an extension or some other composition. 
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Assume that an infinite binary code ziZ2... corresponds to the MSS sequence 
A. Let Ak be the truncation of A so that ziZ2..ifc ^k- One of the "parent" 
branches of ^fe+i is always Aj^. The more distant parent of ^fe+i is denoted by 
yifc+i (and that of A}^ is A}~). For example, the parent branches of A3 = RL^RLR 
are A2 = RL^R and A3 = RL. Ak can be defined also for the case in which Ai. 
lies at the edge of the binary tree — see below. The infinite parents of Afe+i are 
either a{Ak) and e{Ak+i) or e{Ak) and a{Ak+i). It is then clear that A^^i can 
be written either in the form 

or 

Ak+i = ifc+i7fc+i{a(ifc+i)}^'=+' 

where Pk and 7^+1 arc the first binding elements of the proper extension. f3k, hk, 
and nik are needed in the construction of the MSS sequence and can be determined 
by the following rules: 

Rule 1. a) If ik+i — ik, then (3k = {a{Ak)}i+mk and hk is the length of the 
sequence 

a{Ak)A{a{Ak)}^+^^ 
b) mfc+i =mk + hk + l and Ak+i = Ak- 

Proof. With ik+i = ik the same parent branch is approached as on the previous 
step so that Ak+i = Ak- The immediate infinite parent of the new daughter 
branch Ak+i has the form Ak^k^i^k), and the more distant infinite parent has 
the form Akjka{Ak)- Ak+i is their shared beginning. It can be longer than 
Ak = Ak'yk{a{Ak)}T'' only if (3k = {a(ifc)}i+^^. In other words, 

Ak+i = Ak(3k[a{Ak) A {a(Ak)}f+mJ 

lb) can be easily verified by considering the lengths of the sequences in the above 
two equations for Ak-^i, the latter also with the index k. □ 

Rule 2. a) If i^+i ^ i/-, then (3k = {a{Ak-i)} \-\.hk_i and hk is the length of the 
sequence 

a{Ak)A{a{Ak-i)}^^nu-. 
b) m/s+i = hk-i + /ife + 1 and Ak+\ = Ak-\. 
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Proof, ik+i 7^ ^fc entails a turning back towards the k — I'st branch. Thus, 
Ak+i = Af~-i. The infinite parents of the new daughter branch have the forms 
A-k-iPk-i^'iA-k-i) and AkPkO'iA-k)- The shared beginning of the infinite parents 
can be fonger than Ak = ^A;_i/3fe-i{a(^/c-i)}^*'"' only if Pk = {a(^fe-i)}i+/ifc_i- 
Thus, 

Ak+i = AkPk[a{Ak) A {a{Ak-i)}^+h,_,] 

n 

In fact, it would suffice to memorize Ak, hk-i, and ruk because Ak and Ak-i can 
be determined from these. If ii = 0, one has the initial conditions Ai = RLRR, 
ho = 1, and mi = 2 {Ai = R). If ii = 1, then Ai = RLL, ho = 0, and mi = 2 
(ii = 0). 

4. TYPICAL AND EXCEPTIONAL BINARY CODES 

The proofs of Rules 1 and 2 involve only a slight elaboration on the MSS rule. 
The generation of the MSS sequence with these rules is not necessarily much 
more efficient than applying the MSS rule directly. The new formulation helps 
understanding why the length of the MSS sequence increases arithmetically with 
some binary codes and geometrically with others. The increasing length of the 
MSS sequence in a single binary step is given by hk + 1- We call the growth of the 
sequence length arithmetic if sup{Hi, H2, ...} < H < 00, where 

1 ^ 
fc=i 

This definition allows arbitrarily large occasional increments in the MSS sequence 
but they may not be frequent. If the growth is not arithmetic it is called geometric. 

hk must become large for the MSS string length to grow rapidly. Let us denote 
the length of Ak+i by l{Ak+i) and assume that m^+i < /(Afc+i) (call it the simple 
extension condition). Then both Ak and Ak (Rule la) or Ak and Ak-i (Rule 
2a) are long enough so that we can replace the antiharmonic extensions with the 
infinite MSS sequence A in determining hk- This suggests that hk can be large 
only if there is a large block in A identical with the beginning of A and which lies 
after the first 1 + ruk (ik+i = ik) or 1 + hk-i (ik+i 7^ ik) symbols of A. On the 
other hand, the beginning of A can be written as in the form A = Aj^j{a{Aj)}^\.. 
anywhere in the sequence (at arbitrary j). Note that a{Aj) can be replaced by A if 
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l{Aj) > rrij. Thus rrij gives a rough estimate of the length of the above described 
block, identical with the beginning of A, which comes following the first l{Aj) + 1 
symbols. If (ik+i = ik) or hk-i (ik+i ^fc) becomes equal to l{Aj) for some j 
for which rrij is large, then hk can reach a large value. If such an index j is not 
found, then hk can be expected to remain "small". 

The initial values ho and mi are small, and Rule 2b implies that taking steps in 
alternating directions in the binary tree one is not likely to generate large values 
of nij. According to Rule lb, nij grows at least linearly with j if consecutive steps 
are taken in the same direction. In other words, assuming that li — > Aj'yj{A}^^ 
with / Aj we obtain nij > p. In this way we see that arbitrarily large values 
of rrij and, accordingly, of hk are possible. However, this is not sufficient to 
generate geometric growth of an MSS sequence. Geometric growth requires average 
unbounded increases of hk as a function of k. It turns out that judiciously placed 
blocks of identical symbols can bring about such a phenomenon. It is difficult to 
explain the geometric growth of the MSS sequence length for an arbitrary binary 
code. In the following we consider a special form of the binary code whose growth 
properties are easier to understand. 

Proposition. Consider binary codes of the form 1010P^^hOP^^\..10P^''\ Iip{i) > 
and p{i) + < i (z = 2, 3, n), where is calculated from Eqs. (1-5) below, 
then the corresponding MSS sequence has the form 

a-i[ ]-iao[ ]o...an[ ]n 

with an either R or L according to the rule q;„ = cife(n) {k{n) < n) beginning with 
a-i = R and ao = L and 

[ ]n = [ ]fe(n)«j(n)[ ]i(n)«l+j(n) [ ] l+j(n) •••«m(n) [ ]m(n) 

Let l{n) be the length of the block [ ]n and s{n) the length of the whole MSS 

sequence. Then k{n), j{n), m{n), l{n), and s{n) are given by the recursion 
formulae 

s{k{i) - 1) = l{m{i - 1)) (1) 

sij(i) - 1) = lik(i)) (2) 

m{i} ^ j{i}+p{i) -1 (3) 

Id) = simii)) (4) 
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s{i) = s{i - 1) + + 1 (5) 

with the initial conditions k{l) = 0, m(l) = -1, /(-I) = 1(0) = 0, /(I) = 1, 
s{-2) = 0, = 1, s(0) = 2, and = 4. The MSS sequence begins with 

RLL[R]i... 

The proof is a straightforward apphcation of Rules 1 and 2 and is omitted here. 
The condition p{i) + < i is equivalent to the simple extension condition. 

Eqs. (1-4) can be combined into a recursion rule for m{n) alone: 

m(n) = p{n) + m{m{m{n — 1)) + 1) (6) 

In addition to the initial condition for m(l) one has to specify the values 
m(— 1) = m(0) = —2 in order to apply (6). The simple extension condition in 
terms of m{i) becomes m{i) < i — 1. 

With the aid of Proposition binary sequences can be constructed which lead to 
geometric growth of the MSS sequence length. Eqs. (4-5) give 

s{n) = s{n - 1) + s{m{n)) + 1 (7) 

for the length of the sequence. This implies that the window period grows 
according to the recursion g„ = qn-i + 9m(n) with each added block 10^*^"^^ in the 
binary code. The powers p{n) in Eq. (6) can be chosen so that m{n) = n — r with 
r > 1 for n > N. The leading eigenvalue ^ of the transition matrix M defining the 
recursion via (g^, Qn-i, Qn-r+i) = (g^-i, Qn-r) is greater than unity 

and gives the asymptotic growth rate of the MSS sequence length. The simple 
extension condition implies C < 2. 

All Fibonacci sequences correspond to m{n) ~ n — 2 with ( Ki 1.618. Eq. (6) 
implies that the power p{n) takes the constant value 4 after the "transients" have 
died out. 

Example. The MSS sequence for the binary code lOlO^lO^lO^lO^... reads 

RLL[R]iR[RL]2R[RLL{R)]^L[RLL{R)R{RL)]^L[{R)LL[R)R[RL)R{^^ 

Each block a„[ results from a block 10^^"^ in the binary code. The increments 
resulting from h different from zero are in parentheses (...). These increments are 
precisely the same as the blocks [ ]j and m{i) = i — 2 (z = 1, 2, ...). 
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Examples of geometric growth with various r and asymptotic powers p are listed 
in Table I. All were originally found numerically but the third, fourth and fifth are 
beautifully explained by Proposition which predicts correctly the binary period of 
families with the recursion qn = Qn-i +Qm{n) although their binary codes may not 
begin like 101.... More complicated cases can be constructed letting m{n) oscillate 
according to some rule. For example, taking m(2n+l) = 2n—l and m{2n) = 2n— 3 
leads to the recursion relations g2n+i = q2n + q2n-i and g2n = q2n-i + q2n-3 
which can be combined as (j2n+i = ^q2n-i + 92n-3- Prom Eq. (6) one can solve 
for the powers p{2) = 1, p{3) = p{4) = 3, p{5) = 5, p{2n) = 4 (n = 3,4,...), 
and p{2n + 1) = 6 (n = 3,4,...). Each asymptotic binary block 10^10^ means 
multiplying the period with the average factor C ^ 1.554. Aperiodic binary codes 
can be generated making the oscillations in m(n) aperiodic. Instead of letting m(n) 
oscillate regularly like above, we took m(n) = 2[n/2] — 3 or m{n) = 2[n/2] — 1 
at random ([ ] stands for the integer part), beginning with an adjustable n. The 
longer the leading regular part, the faster the MSS sequences grew in length. On 
account of inherent limitations of computers, no certain conclusion could be drawn 
about the asymptotic algebraic or geometric growth of such sequences. It may be 
interesting to notice that the sum of the second differences of the MSS lengths 
seldom seemed to display a sustained growth suggesting an asymptotic algebraic 
fate to all MSS sequences arising from random period binary codes. 

A purely numerical method of finding geometric growth consists of first defining 
a recursive rule for the increase of the window period and then determining the 
binary code which yields the windows (only a small number of rules lead to simple 
behavior). Each "appropriate" recursion rule has a characteristic repeating pattern 
of the periodic binary tail. For example the rule qn = q-n-i + qn-2 gives the same 
repeating pattern 10^ for all choices of qi and q2 (if one always moves to the 
right in the binary tree). Taking off from an arbitrary window, however, and 
successively adding the pattern 10^ to the binary code will almost certainly lead 
to asymptotically arithmetic increase of the MSS sequence length! 

Geometric growth obviously requires a very synchronous binary code. A single 
mismatched binary symbol may suffice to turn the growth into arithmetic. 
Therefore, one expects the MSS sequence length to grow arithmetically for most 
binary codes with periodic tails. We call such codes "typical" whereas the codes 
leading to the geometric growth of the MSS sequence are called "exceptional". 
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The examples of Table I have been constructed by always picking the closest 
matching window on the right hand side. We could not find any exceptional codes 
with periodic tails by choosing the next window on the left hand side. One may 
begin, nevertheless, with a binary code that carries a number of leading zeros. 

5. NUMERICAL SCALING RESULTS 

The scaling behavior is strongly related to the manner of growth of the length of 
the MSS sequence along a path in the binary tree. In the following, the scaling 
of both typical and exceptional binary codes is discussed concentrating on codes 
with periodic tails. 

5.1 Typical codes with periodic tails 

The binary sequence 1°° corresponds to the stable periods 4, 5, 6, ... approaching 
the final crisis point. The scaling properties of this sequence are well understood 
[10,11]. Along the parameter axis the scaling of the superstable parameter values 
is geometric. The scaling factor is determined by the derivative of the map at the 
unstable fixed point which is the image of the critical point under the second iterate 
of the map. The widths of the windows scale by the square of this factor. Every 
other binary code with the tail 1°° corresponds to a window sequence approaching 
a tangent bifurcation point, i.e. the left end point of some periodic window. It 
has been shown [12] that this results in slower than geometrical scaling. 

On the other hand, every binary sequence with the tail 0°° leads to a sequence 
which accumulates at an internal crisis point. All these points are fully chaotic 
according to a theorem by Misiurewicz [1,4]. In the same way as for the final 
crisis point, the scaling properties are determined by the Lyapunov factors of the 
associated unstable orbits (Section 6). 

Let us now consider an arbitrary binary code of the form /J°° and let /x^ and 
A/z/c be the superstable parameter value of the window IJ^ and its width. We 
determine the scaling factors 

/ifc - ^^k-l A//fc-i 

~ j (^k — — 7 

fJ'k+l — IJ'k i-^l^k 

for the logistic map and find the period v by which 5k and oscillate as A; — > oo. 
The asymptotic limits of the products of v subsequent 5k ^ and u/s's are denoted by 
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6 and a, respectively. Table II displays f and 6 for a number of examples. Neither 
v nor S is universally determined by the tail of the binary code. In Section 6 it is 
shown that both are related to the orbit of the critical point at the accumulation 
point of the window sequence. The accumulation point for this class of codes is 
always a Misiurewicz point at which the critical point is mapped onto an unstable 
period Mp after Mi iterations. A proof of this claim is given in Appendix. The 
numbers Mp and Mi are included in Table II. It turns out that an MSS sequence 
corresponding to the typical binary code has a periodic tail, which is possible only 
if the critical point is mapped onto an unstable orbit. The repeating pattern gives 
the kneading sequence of the unstable orbit. The quantity u is greater than unity 
if the whole pattern is not traversed in one period of the binary tail. It may happen 
that the pattern is run through more than once per binary period. In Table II, p 
gives the number of times the pattern is completed in ly periods of the binary tail. 

In Section 6 the result of Post and Capel [11] that u = 5^ is generalized to 
all codes within this class. In other words, the scaling behavior in this class is 
equivalent to the one for the code 1°° and those with the tail 0°°. 

5.2 Exceptional codes with periodic tails 

It has been shown in the Fibonacci case that the scaling of the superstable 
parameter values is superexponential [13,14]. A geometric increase of the MSS 
sequence length leads to superexponential scaling in all the cases we studied. 
Furthermore, a comparison of the appropriate scaling factors [13,14] for the logistic 
map on the one hand and f^{x) — //sin(7ra;) on the other suggests that the 
superexponential scaling is universal. 

5.3 Codes with aperiodic tails 

It was demonstrated in Section 4 that exceptional aperiodic codes with rapid 
growth rates of the MSS sequence length can be explicitly constructed. For 
a random code with no build-in blocks of identical symbols, however, it is 
natural to expect arithmetic average growth. The numerical studies bear out this 
expectation. We observe average geometric scalings of the superstable parameter 
values and positive Lyapunov exponents at the accumulation points. 
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6. GEOMETRIC SCALING 



We are interested in families of MSS sequences whose binary codes have periodic 
tails. These sequences have some initial length Mj. They grow in steps of n = pMp 
per advancing v periodic blocks of the binary tail. In order to get at the scaling 
parameters 5 and a defined in the preceding section, we take advantage of the 
local description of Post and Capel [11]. It describes the map fjj^{x) within the 
window of period I. We choose a particular representation of the logistic map 
ffj-i^) = 1^ ~ where z may take the values 1 < ^ < oo (i.e., in this section 
wc do not restrict ourselves to the case of a quadratic maximum). Then the local 
descriptions is centered at x = 0, and the following reduced form applies 

= Piil^) + Xilxiil"" + Ki\xii\'^^ + ... (8) 

where a;/(j=o) is some starting point close to the central value x = 0, so close that 
we do not have to worry about the term to the power 2z on the right hand side. 
The local description then incorporates one parameter A;. Finally 

piifi) = ii-\ii-\ii- (9) 

with / — 1 pairs of vertical bars. For the present purposes the index / will need 
to take values like m + kn and m + {k + l)n where A; is a large integer since we 
are interested in the asymptotic /x scaling and the ratios of widths of consecutive 
windows in our families. 

With the substitutions 

xii — Ui\X\~^^^^~^'^ sgnX 

p= -r\X\-^/^^-^hgnX (10) 

and leaving out the subscript of A, the above form turns into Post's and Capel's 
normalized submap 

Ui+i = |tti|^ - r (11) 



First consider the ratios of window widths. If p varies by Ap when p, varies across 
the period-/ window, the window width expressed in p, is obviously 

dp 



dp, ' 
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where a tilde on the symbol /i indicates the superstable parameter value. By the 
scaling, Eq. (10), between p and r we can express Ap in terms of Ar, the invariable 
normalized window width of Eq. (11). The "physical" window width wi is then 

Ar 

wi = Ap 



A|i/(^-i)(4^| ) 



We now give the index I the values N + n and N where N = m + kn. For ratios of 
window widths, Ar cancels, and the scaling in the family becomes asymptotically 



X 1/{Z-1)^\ 

= WN/WN+n = ) — (12) 



when = m+nk becomes large, i.e. k becomes large (at is essentially the product 
of v subsequent a^s defined in Section 5). 

For the scaling of the positions of the windows we turn to the function Pm+kn- It 
determines the height of the quadratic maximum or minimum of the local mapping 
jm+kn^rj.'^ and it gives the value of x to which the center point a; = is sent in 
jm+kn^rj.y At p = Pm+nk, Pm+nk = p(a*) takes the valuc zero. Let us look at 
the functions Pm+nk for different k at the accumulation point poo of the family. 
Denote again N = m + nk. From the definition of p above, Eq. (9) 

Piv(/^) = 4^(o) 

and 

But the orbit of the center x = repeats itself with the period n. It follows 

PN+n{Poo) - Pn{Poo) 

which means that the Pm+kn for different k all meet at the same point at poo- 
Assuming that the p{p) are locally straight lines, we get the asymptotic p scaling 
from the derivatives of the p for different k. Estimate pn and pN+n as follows 

Pn{Poo)^—^\ {pN - Poo ) 

and 

/ N dpN+n I / ~ X 

PN+n{P^) ^|p=/i^+„ W+n - P<^) 
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These two quantities being equal, 



dpN+n I 



jiN+n - f^oo ^ I - 



and 



/dpN+n\ \i/dpN\ \ (.r,\ 

again for large N. We have not seen this result in the literature. 

What remains is calculating the scaling factors A and the ratio of {dpN+n/dfJ') at 
fi = p,N+n and dppf/dfj, ai jj, = p,N . We follow Post and Capel [11]. 

Take A first. In the following, a tilde on the variable x refers to an iterate of the 
central value x = 0. At superstability, since p vanishes there, we have from the 
local description of Eq. (8) 

xi(^i+i^ = X\xii\'' (14) 

where some starting point xu close to the central value x = has been picked. 
Wc will express the left hand side of this equation as an expansion in terms of the 
following quantity 

V = ffli {Xii)- fjXi (xo) 

with the same point xu . In a single shot of the mapping /^j, the central point 
and a point in its vicinity will be sent far away from the center but roughly to the 
same location. Therefore the quantity v is small. To linear order in v then 



= /t'(/M,(^o)) + ff ^' {U^{xo))v = n (15) 

i=i 

where a prime denotes a derivative with respect to the argument of the function 
at the indicated value of the argument. For the specific form f{x) = p — \x\^ 



V = -\xu 

and Eq. (15) yields 



I 2 



l-l 

xiii+1) = -Ylfp,/{xj)\xii\ 
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and one reads from Eq. (14) 

l-i 

^=-\[h/(5:,) (16) 



With cleared, one still needs to calculate the derivative of p with respect to ji 
in Eqs. (12) and (13). Look again at the reduced map /^(x) 

xi = p + A|a;o|^ 

and its derivative with respect to jj, 

dxi dp ^ ^1 1^ 
dji dp, dp, 

Choose xq as xo = and take the expression at the superstable point. The second 
term on the right vanishes and 

dxi I dp I 
dp, ~ dp, 

One thus needs an expression for (dxi/dp,)-^. Consider the map 

In general depends on xq in addition to p which is indicated by the partial 
derivatives in the sequel 



d^i-'rl -, I 12 — 1 d^i 1 , /• // \dXi 

— ^ — = 1 — z\Xi\ sqnxi—— = 1 + /"„ (xj— — 
dpL ^ ^ ' dpi ""^ ^ ' dpi 



= 1 + f^'{xi)[\ + //(x,_i)[l + //(a;,_2)[l + • • • [1 + f^ixM ■ ■ •] (17) 

Nothing prevents us from choosing i+l=l and picking the superstable pii along 
with = = 0. Then the iterates depend only on pi and 

^1^, = 1 + /m/(^^-i)[1 + /a/(^/-2)[1 + - • • [1 + /a/(^i)]] • • •] = 1 + E n 

»=1 j=i 



or finally 

dp. 



Mi 
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Now return to Eq. (12). From Eq. (16) we get immediately at the limit when 
N ^oo 

K+n 

AiV+n/AjV = n ^'-^N+S^i) 
i=K 

where K is some number larger than Mj, the number of initial iterations before 
hitting the unstable period, and smaller than N — Mf, Mf being the number of 
the last steps drifting away from the unstable period to end at the center. It is 
important to notice for what follows that K may vary over the whole middle range 
of the N + n cycle. It turns out convenient to call the derivatives in the product 
and L the whole product from i = K to i = K + n. The contribution to ak of the 
first factor in Eq. (12) is then 

\ 1/(2-1) 

^Xn+u^ = i:i/(-i) (18) 

The second quantity in Eq. (12), 

f^l )/A ) (19) 

has a structure which is easiest grasped looking at one derivative at a time. In the 
following the products are written in the order of increasing length, more or less 
like in Eq. (17). We need clarifying notation: 

1 + Ln-i + Ln-iLn-2 + Ljv-iivjv-2-^Ar-3 + ■ ■ ■ + e — E 

where e = Ln-iLn-2Ln-z ■ • -L^-Mf Then 
dp I 

-r ,,-r, = 1 + Ln-1 + + Ln-iLn-2Ln-s + H Ln-iLn-2 • • ■ Li 

= E -\- e(Ljv_M/-l + LN-Mf-lLn-Mf-2 + ■ ■ ■ 
+LN-Mf-lLN-Mf-2 ■ ■ • LM-Mf-n+1 + L) 
+eL{LN-Mf-n-l + LN-Mf-n-lLN-Mf-n-2 + ' ' ' 
+LN-Mf-n-lLN-Mf-n-2 ' ' ' LN-Mf-2n+l + L) -\ (20) 

Now the expressions in the brackets in Eq. (20) are identical since they run over 
full unstable periods of length n. We again introduce a new symbol B, this time 
for the sum in the brackets multiplied by e. Then 

^1 . =E + B + BL + BLL + --- + BL^''-^f-'^'y''+P (21) 

dlJ,^tJ-=fJ-N 
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where 



P = eL(^-^/-^^)/-(LM,-i + Lm,-iLm,-2 + ■■■ + Lm,-i ■ ■ ■ Li) 

The quantity Mf may be thought of as having been chosen such that {N — Mf — 
Mi)/n is an integer. P represents the in-going steps before hitting the unstable 
period. 

Introduce yet another symbol S = S{N) for the right hand side of Eq. (21) 
without the term E. Remember that the ratio we are seeking to calculate, Eq. 
(19), is between two expressions of type Eq. (21) with the numerator having N + n 
in the place of the denominator's N. If the denominator is expressed as in Eq. 
(21), the numerator cycle is longer by n and the right hand side of Eq. (21) will 
have one more term, 

j^j^{N+n-Mf-Mi)/n 

and the factor P becomes multiplied with L. This can be interpreted as multiplying 
by L the terms which were given the name S above and adding back a B. It follows 

S is obviously rapidly growing function of provided that L is larger than unity, 
i.e. that the periodic orbit at the accumulation point of the family is unstable. 
All other quantities in Eq. (22) remain constant when N grows. The asymptotic 
result at the limit A?" — > oo is simply 

.dp, y.dp, , ^ 

^dp ^l^=P'N+n'' \n=iiNl 

By Eq. (13) this result is directly the parameter scaling factor 5. 



Together with Eq. (18) we get for the scaling of the window width, Eq. (12) 



7. DISCUSSION 

It was shown in Section 6 that typical binary codes with periodic tails lead to 
geometric scaling whatever the order z of the critical point. The situation is very 
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different in the case of the exceptional codes. Our numerical results suggest (see 
Section 5) that the scaling for an exceptional code with periodic tail is always 
superexponential in the case of the logistic map with a quadratic maximum. If we 
let the order of the maximum increase, however, the scaling becomes geometric and 
universal. For the Fibonacci sequences this has been pointed out by M. Lyubich 
[19]. The critical value of the exponent z for the Fibonacci sequences is 2. In 
other cases it is numerically much harder to pinpoint the exact critical value of 
the universality class at which the transition from superexponential to geometric 
scaling takes place. This would be an interesting problem for future studies. 

A binary tree of stable periodic attractors appears also in invertible circle maps 
[20]. It is therefore interesting to compare the scaling behaviors in these two 
cases. In a circle map the period of the attractor increases geometrically for almost 
all routes in the binary Farey tree. In the standard case this leads to universal 
geometric scaling. In the unimodal map a geometric increase of the period and 
the resulting universal scaling is observed only in some exceptional cases. It is an 
open question whether one could construct a "non-renormalizable" circle map with 
superexponential scaling. It is not clear how the various scaling properties relate 
to the Lebesgue measure of the parameter values for aperiodic attractors. For a 
standard critical cubic circle map this measure is zero whereas for the quadratic 
unimodal map the measure of aperiodic attractors is positive. 

The typical binary codes with periodic tails give a huge number of new 
Misiurewicz points in a unimodal map (recall that the same structure is repeated 
within each window). These points form a subset of the set of all Misiurewicz 
points of the Mandelbrot set [21]. It would be an interesting problem in the 
theory of numbers to work out whether the typical binary codes with periodic 
tails (taking into account the self-similarity) give all the Misiurewicz points of a 
unimodal map. 

APPENDIX. A TYPICAL BINARY CODE WITH A PERIODIC 
TAIL LEADS TO A MISIUREWICZ POINT 

Let us assume that the binary code has a periodic tail which is neither 0°° nor 
1°°. The code can always be written in the form IJJJ...^ where / and J = ji-.-jx 
are finite binary codes with j2 = ji. Let A be the corresponding infinite MSS 
sequence. 
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Lemma. Consider a binary code of the above form with < H < oo for every k. 
Then the accumulation point of the corresponding sequence of periodic windows 
is a Misiurewicz point. 

Proof. Because j2 = ji, Rule 2b implies that mN+2+qK < 2-ff for q = 0, 1, 2, 
where N is the length of /. Furthermore, by Rules lb and 2b rrik+i < Tnk + H. 
Combining these two results gives an upper bound for mk- irtk < [K + 1)H for 
k — N + 2, + .... Thus, there exists a,p > N + 2 such that nip < l{Ap). Because 
mk+i < l{Ak+i) for k = p— 1, p, it is not necessary to carry out the antiharmonic 
extensions when determining hk {k > p — 1) in Rules la or 2a. In particular, A can 
be written as the composition Ap-iPp-i{A}i^^^ Pp{A}i^ .... The set of possible 
values for nik and hk is finite when k > N + 1. Therefore, there exist finite P > p 
and Q so that {mp,hp,hp-i) = {mp+QK,hp+QK,hp+QK-i). Because ip+k = 
ip+QK+k for k = 0,1,2..., Rules 1 and 2 imply that {mp+k, hp+k, hp+k-i) = 
{mp+QK+k, hp+QK+k, hp+QK+k-i) which is possible only if A has a periodic tail. 
By Corollary II. 8. 4 of ref. [1], A is an MSS sequence for a Misiurewicz point. □ 
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Table I. The repeating patterns of the periodic binary tails for some recursive 
rules leading to an asymptotic geometric increase (given by the factor Q of the 
MSS sequence length. 



rule 


c 


pattern 


Qn — '^Qn-l + Qn-2 


2.414 


10^^ 


Qn = Qn-1 + '^Qn-2 


2.000 


10^ 


Qn = Qn-1 + Qn-2 


1.618 


10^ 


Qn = Qn-1 + Qn-3 


1.466 


10^ 


Qn = Qn-1 + Qn-4: 


1.380 


108 


Qn = Qn-2 + Qn-3 


1.325 


10^0 



Table II. The scaling period ly and the factor 6 for some typical binary codes with 
periodic tails. Mp gives the period of the unstable orbit and Mi the length of the 
transient at the accumulation point. The repeating MSS pattern is gone through 
p times during u periods of the binary tail. 



code 


u 


P 


6 


Mp 


M, 


(01)°° 


1 


2 


6.996 


2 


4 


(10)°° 


1 


2 


3.716 


1 


4 


0(01)°° 


1 


2 


5.560 


2 


6 


1(10)°° 


1 


2 


3.931 


1 


5 


(001)°° 


1 


1 


12.11 


6 


4 


(010)°° 


2 


1 


53.96 


9 


10 


(100)°° 


1 


1 


4.962 


3 


8 


(110)°° 


1 


1 


7.694 


3 


4 


(101)°° 


1 


1 


6.736 


3 


3 


(011)°° 


2 


1 


95.30 


9 


4 
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FIGURE CAPTION 

Fig. 1. The beginning of the infinite binary tree of periodic windows. 
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